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ABSTRACT 

We present results from strong-lens modelling of 10,000 SDSS clusters, to estab- 
lish the universal distribution of Einstein radii. Detailed lensing analyses have shown 
that the inner mass distribution of clusters can be accurately modelled by assuming 
light traces mass, successfully uncovering large numbers of multiple- images. Approx- 
imate critical curves and the effective Einstein radius of each cluster can therefore 
be readily calculated, from the distribution of member galaxies and scaled by their 
luminosities. We use a subsample of 10 well-studied clusters covered by both SDSS 
and HST to calibrate and test this method, and show that an accurate determina- 
tion of the Einstein radius and mass can be achieved by this approach "blindly", in 
an automated way, and without requiring multiple images as input. We present the 
results of the first 10,000 clusters analysed in the range 0.1 < z < 0.55, and compare 
them to theoretical expectations. We find that for this all-sky representative sample 
the Einstein radius distribution is log- normal in shape, with {LogiOe")) = 0.73to.o3' 
(J = 0.316lo"oo2 5 with higher abundance of large Oe clusters than predicted by 
ACDM. We visually inspect each of the clusters with 0^ > 40" {zg = 2) and find that 
~ 20% are boosted by various projection effects detailed here, remaining with ~ 40 
real giant-lens candidates, with a maximum of 0^ = 69 ± 12^' {zg = 2) for the most 
massive candidate, in agreement with semi- analytic calculations. The results of this 
work should be verified further when an extended calibration sample is available. 

Key words: cosmology: theory, dark matter, galaxies: clusters: general, galaxies: 
high-redshift, gravitational lensing: strong, mass function 



1 INTRODUCTION 

Clusters of galaxies play a fundamental role in testing cos- 
mological models, by virtue of their position at the high 
end of the cosmic mass spectrum. Massive galaxy clusters 
gravitationally-lens background objects, forming distorted, 

^ E-mail:adiz@wise. tau.ac.il 



magnified, and often multiple images of the same source, 
when the cluster surface density is high enough. These ef- 
fects are in turn used to map the gravitational potentials 
and mass of the lensing clusters, hence providing some of 
the best constraints on the nature and shape of the under- 
lying matter distributions (Broadhurst et al. 2005a, Bradac 
et al. 2006, Coe et al. 2010, Zitrin et al. 2010, Merten et al. 
2011). 
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Large sky surveys such as the Sloan Digital Sky Survey 
(SDSS; see Abazajian et al. 2003,2009) ahow for important 
scientific work with different astrophysical implications (e.g., 
Tegmark et al. 2004, 2006, Tremonti et al. 2004, Eisenstein 
et al. 2005, Seljak et al. 2005, Wojtak, Hansen, & Hjorth 
2011). The large amount of data enables extensive studies 
with a clear statistical advantage. Here we make use of the 
results of a new cluster-finding algorithm operated on the 
SDSS DR7 data (Hao et al. 2010; on DR7 data see Abaza- 
jian et al. 2009), in order to derive the Einstein radius dis- 
tribution of a significant, statistical sample. As presented in 
their work, more than 55,000 clusters were found using this 
successful and rather conservative algorithm, which we have 
taken upon to analyse using our improved lensing-analysis 
tools (e.g., Zitrin et al. 2009b, see more details in §2), pre- 
senting here the results of the first 10,000 clusters analysed. 

The effective Einstein radius plays an important role 
in various studies. The Einstein radius describes the area 
in which multiply-lensed images may be seen due to the 
high mass-density of the cluster. By definition, within this 
critical area the average mass density is equal to T^crit (for 
symmetric lenses), the critical density required for strong- 
lensing, whose value is dependent on the source and lens dis- 
tances. In general, obtaining the critical curves with great 
accuracy allows matching up multiple-images, which in turn 
help to improve and better-constrain the model in order to 
derive the mass distribution and profile more accurately, 
teaching us about certain properties of both the observed 
and unseen matter. The Einstein radius therefore consti- 
tutes a measure of the strong- lens size (and efficiency), and 
directly enables us to estimate the amount of mass enclosed 
within it; Be = (^ 4GM(^<0e) ^^^1/2 symmetric lenses (e.g., 
Narayan & Bartelmann 1996, Bartelmann 2010), where di^ 
ds, and dis, are the lens, source and lens-to-source (angular- 
diameter) distances, respectively. Equivalent ly, the effective 
Einstein radius used here is simply a measure of the critical 
area. A, so that 6e = \J Ajix. 

In recent years it has been proposed that the Einstein 
radius distributions of several small samples of clusters, pose 
a challenge to ACDM (e.g., Broadhurst h Barkana 2008, 
Zitrin et al. 2009a, 2011a). Other discrepancies such as the 
arc abundance, several uniquely large Einstein radii, mas- 
sive high-z clusters, high NEW concentration parameters, 
and comparison to N-body simulations, contribute further 
to this tension, though most studies show mainly a mod- 
erate discrepancy (e.g., Bartelmann et al. 1995, Wambs- 
ganss et al. 1995, Dalai, Holder, & Hennawi 2004, Broad- 
hurst et al. 2005b, 2008, Hennawi et al. 2007a,b, Hilbert et 
al. 2007, Sadeh h Rephaeh 2008, Oguri h Blandford 2009, 
Oguri et al. 2009, Puchwein & Hilbert 2009, Meneghetti et 
al. 2010a,2011, Sereno, Jetzer & Lubini 2010, Gralla et al. 
2011, Horesh et al. 2011, Umetsu et al. 2011a, Zitrin et al. 
201 la, c). Obtaining a credible empirical distribution of Ein- 
stein radii from an unprecedentedly large sample is of clear 
value, welcoming in addition complementary mass measure- 
ments through similarly automatic weak-lensing analyses 
(e.g., Hildebrandt et al. 2011) and other observations, such 
as of X-ray emission or the SZ-effect, when possible. 

The advances in computational power over the past 
decades along with higher quality data and our efficient 
method for analysing strong-lenses (Broadhurst et al. 2005a, 
Zitrin et al. 2009b) now enable such an extensive study. 



Based on previous analyses of many clusters, we now se- 
curely determine typical physical parameters to which the 
critical curves are relatively indifferent, so that we extrapo- 
late and test these assumptions to perform our analysis on 
the sample presented here. In particular, in this work we de- 
scribe a simple and efficient method to model cluster-lenses 
based on the light distribution of bright cluster members, 
which as we have targeted to show, allows to derive the Ein- 
stein radius with sufficient accuracy, in an automated mode. 

Automated surveys for lensing have been presented be- 
fore, though mostly based on the observed arc properties, 
or relate to either galaxy- lensing scale or the weak-lensing 
regime (e.g., Webster, Hewett h Irwin 1988, Cabanac et al. 
2006, Mandelbaum et al. 2006, Johnston et al. 2007b, Cor- 
less & King 2009, Marshah et al. 2009, Sheldon et al. 2009, 
Bayliss et al. 2011a,b, Hildebrandt et al. 2011), and have yet 
to produce statistically-significant results for the Einstein 
radius distribution directly from SL modelling. Other avail- 
able SL methods, though can be successful, either require 
the location of many multiple-images as input or currently 
have too many free parameters, rendering such a "blind" 
study impossible. 

The SL modelling method we implement here is based 
on the reasonable assumption that light approximately 
traces mass, which we have shown is most efficient for finding 
new multiple-images as the mass model is initially well con- 
strained with sufficient resolution to derive well-approximate 
critical curves (see Broadhurst et al. 2005a, Zitrin et al. 
2009b, 2011a,b,c, Merten et al. 2011). Recently we have 
tested the assumptions of this approach in Abell 1703 (Zitrin 
et al. 2010), by applying the non-parametric technique of 
Liesenborgs et al. (2006, 2007, 2009) for comparison, yield- 
ing similar results with only minor differences in the over- 
all mass distribution and critical curves, especially where 
galaxies are seen since they are not included in the non- 
parametric technique. Independently, it has been found that 
SL methods based on parametric modelling, i.e., based on 
physical assumptions or parametrisations (for other para- 
metric methods see, e.g., Keeton 2001, Kneib et al. 1996, 
Gavazzi et al. 2003, Bradac et al. 2005, Julio et al. 2007, 
Halkola et al. 2008), are accurate at the level of a few per- 
cent in determining the projected inner mass (Meneghetti et 
al. 2010b). Clearly, non-parametric techniques and methods 
that are based directly on arc morphologies are also impor- 
tant: non-parametric techniques (e.g., Diego et al. 2005, Coe 
et al. 2008, Merten et al. 2009) are novel in the sense that 
they are assumption- free and highly flexible (e.g., Coe et 
al. 2010, Ponente & Diego 2011), and methods based di- 
rectly on arc morphologies yield high resolution results (see 
also Grillo et al. 2009). The parametric method presented 
here, is simply aimed to produce the critical curves in an au- 
tomated way based on simple physical considerations (and 
thus is capable of flnding multiple images as we have shown 
constantly before), and constitutes another important step 
towards the ability to deduce the lensing properties of clus- 
ters in large sky surveys in an automated way, so that we aim 
now to present the first observationally-deduced, universal 
distribution of Einstein radii. 

The incorporated method involves only four free pa- 
rameters. Three of them are known sufficiently well a-priori 
and have only negligible effect on the critical curves and 
resulting Einstein radius, for which we adopt typical val- 
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Figure 1. Calibration sample. Einstein radii (for Zs = 2) derived 
by our "blind" automated algorithm in SDSS data and based on 
the assumption that light traces mass, with a typical M/L, versus 
the Einstein radii derived by detailed analyses of HST images of 
the same clusters and using the multiple images as constraints 
(e.g., Broadhurst et al. 2005a, Richard et al. 2010, Zitrin et al. 
2010, 2011a,b). As can be seen, the "blind" method, based on the 
light distribution of bright cluster galaxies and without using any 
information regarding the location of multiple-images, shows re- 
markably similar results to those derived by the detailed indepen- 
dent analyses (see also Fig. |14| ). The errors in the Einstein radii 
are typically ~ 10%, and overplotted is also an x = y dashed line. 
The comparison sample spans the redshift range 0.15 < zi < 0.55. 
Blue open circles are the results from the blind analysis with the 
best-fitting parameters derived from minimising by all 10 clusters 
together, while red filled circles are the results from a " Jackknife" 
minimisation described in §3.2.2| in order to demonstrate how 
well one could assess the Einstein radius, we perform the min- 
imisation for 9 different clusters at a time and analyse the tenth 
cluster with the resulting parameters. With this we obtain devia- 
tions of up to ~ 17% in our ability to blindly estimate the critical 
curves by the automated procedure described in this work. In a 
complementary error-propagation check ( j3.2.2| > we obtain simi- 
lar results of Icr ~ 18%, which we take hereafter as the errors for 
the full-sample analysis. 



ues deduced from detailed analyses of a few dozen clusters 
(more details are given in The fourth parameter, which 
varies from cluster to cluster, is the overall (mass) normal- 
isation, but since the respective distances are known, this 
can be simply overcome by finding a typical mass-to-light 
ratio (M/L) normalisation. The M/L term is embedded, 
in practice, in a redshift-dependent normalisation factor, 
which is iterated for the best fit using 10 clusters which have 
been accurately-analysed in HST images and have parallel 
SDSS data listed in the Hao et al. (2010) catalog. These 
include some well-known lensing clusters such as A 1689, 
A1703, MS1358, Z2701, and others (see, e.g., Broadhurst et 
al. 2005a, Richard et al. 2010, Zitrin et al. 2010, 2011a,b). 
The results of this comparison are shown in Figure [l] and 
Table [U 

The paper is organised as follows: In §2 we detail the 
modelling and the assumptions on which our algorithm is 
based. In §3 we discuss the results and relevant uncertainties, 
which are then summarised in §4. Throughout this paper we 
adopt a concordance ACDM cosmology with (Qmo = 0.3, 




Figure 2. The general starting point of our lens model, where we 
define the surface mass distribution based on the cluster member 
galaxies (see Q listed in the Hao et al. (2010) cluster catalog. 
In this figure we show the lumpy (galaxy) component for Abell 
1703 as an arbitrary example (see also Zitrin et al. 2010 for an 
equivalent figure but from HST observations. Axes are in pixels 
with 0.2'' /pixel). We perform the same simple procedure for each 
of the 10,000 clusters drawn from the Hao et al. catalog. 



Qj^Q — 0.7, h — 0.7). All Einstein radii referred to in this 
work are for a fiducial source redshift oi Zs = 2. We also 
note that all logarithmic quantities in this work are in base 
10, unless stated otherwise, and are denoted conventionally 
as ''Log". 



2 STRONG-LENS MODELLING AND 
ANALYSIS 

The method we apply here is based on the simple assump- 
tion that mass traces light. This well-tested approach to 
lens modelling has previously uncovered large numbers of 
multiply-lensed galaxies in ACS images of e.g., Abell 1689, 
C10024, 12 high-z MACS clusters, MS1358, "Pandora's clus- 
ter" Abell 2744, and Abell 383 (respectively, Broadhurst et 
al. 2005a, Zitrin et al. 2009b, 2011a,b, Merten et al. 2011, 
Zitrin et al. 2011c). As the basic assumption adopted is that 
light approximately traces mass, the photometry of the red 
cluster member galaxies is used as the starting point for the 
mass model. 



2.1 Initial Mass Distribution 

We now wish to calculate the deflection field by the cluster 
galaxies, or the initial mass distribution. By assuming that 
the flux is proportional to the mass, i.e., assigning a certain 
M/L ratio, the deflection field contributed by each galaxy 
can now be calculated by assigning a surface-density profile 
for each galaxy, S(r) = Kr~^, which is integrated to give 
the interior mass, M(< 0) = ^^{diO)'^~'^ . This results in a 
deflection angle of (due to a single galaxy): 



a{0) 



4GM(< 



c^O dsdi ' 



(1) 



or more explicitly by inserting M(< 6) from above: 
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Figure 3. Smoothed mass distribution. To represent the DM dis- 
tribution we smooth the lumpy component (Fig.[2j of each of the 
10,000 clusters drawn from Hao et al. clusters catalog (see also 
Zitrin et al. 2010 for an equivalent figure based on HST observa- 
tions. Axes are in pixels with ^.2" /pixel). This smoothing proce- 
dure is most useful in generating, when combined with the lumpy 
component, a very reliable deflection fleld and corresponding crit- 
ical curves, as we have shown for many clusters (e.g., Broadhurst 
et al. 2005a, Zitrin et al. 2009b, 2011a,b,c, Merten et al. 2011), 
allowing to flnd large numbers of multiple-images by the model. 



a{9) = 



2 — q L 



dis 



(2) 



We note that all quantities are known, except for 
the normalisation factor which is related to the M/L ratio 
(note that q is maintained constant on a typical and known 
value, see {2.3). Thus, finding the explicit term for K which 



scales correctly all clusters we analysed to date (taking into 
account the different lens and source redshifts) allows us - 
in principle - to perform the automated survey of Einstein 
radii, following the procedure described below. 

By defining Kq — ^di'^~^ we can reduce the 

latter formula to get: 



a(e) = Kq 



Ql-q 



(3) 



where Kq also depends on the redshifts involved, and on the 
power- law index, q (which is set to constant throughout. 

The deflection angle at a certain point due to the 
lumpy galaxy components is simply a linear superposition 
of all galaxy contributions scaled by their luminosities, Li 
(in L0 units): 



agal{0)=KqJ2T^\0-0^\ 



\e-e^ 



(4) 



In practice we use a discretised version of equation |4] 
over a 2D square grid Om oi N x N pixels, given by: 



O^gal,x{0m) — Kq 

i 



Li Atjmi 



(5) 
(6) 



where (Axmi, A^/mi) is the displacement vector Om — Oi of 
the mth pixel point, with respect to the ith galaxy position 

Note that to obtain the luminosity Li of each member, 
we convert its SDSS r-band luminosity to the corresponding 
(Vega) B-band luminosity by the LRG template given in 
Benitez et al. (2009). 

From these expressions a deflection field for the galaxy 
contribution is easily calculated analytically as above, and 
the mass distribution is now rapidly calculated locally from 
the divergence of the deflection field, i.e., the 2D equivalent 
of Poisson's equation. An example is given in Figure [2] 

2.2 The Dark Matter Distribution 

The mass contribution of galaxies is anticipated to comprise 
only a small fraction of the total mass of the cluster, which 
is expected to be dominated by a smooth distribution of 
DM. We now simply assume that the galaxies approximately 
trace the DM. As mentioned, this assumption was found to 
work very well in earlier work on many clusters where large 
numbers of multiple-images were found accordingly. These 
multiply-lensed systems are not simply eye-ball candidates, 
but are reproduced and predicted by the preliminary model, 
indicating that this model, based on the assumption that 
light traces mass, is initially well constrained. 

Since the DM is of course expected to be smoother than 
the distribution of galaxies, we smooth the initial guess of 
the galaxy distribution obtained above, choosing for conve- 
nience a low-order cubic spline interpolation, typical to the 
many previous analyses mentioned above. The smoothing 
degree (the polynomial degree, S) is also a free parameter 
of the model, and the deflection field contributed by the DM 
is then simply the sum of the contribution from each point 
(or pixel) in this smooth DM component. This smoothing 
procedure is the key to our method's success in locating 
multiple-images, and is in practice more useful than assum- 
ing a general DM shape such as NFW or pseudo-isothermal 
spheres, which are highly symmetric and do not necessar- 
ily describe the complex inner DM distribution in detail, 
often not allowing to find in advance the multiple images 
according to the initial mass distribution. An example of a 
smoothed component is shown in Figure [S] 

The deflection field of the DM is then (where each pixel 
is treated as a point mass) given by: 



AXn 



[{Ax^^)^ + {Ay rr 



CUDM.yiOm) = Kg ^ Pi 



Ay„ 



[{Ax^,)^ + [Ay^iy] ' 



(7) 



(8) 



where Pi represents the (unnormalised) mass value in the ith 
pixel of the smooth component. We therefore obtain now the 
deflection field due to the DM, hereafter aoMiO). 



2.3 The Total Deflection Field 

Having calculated the two components of the deflection field, 
we now simply combine them to get a total deflection field 
as follows: 



ariO) = Kgaidgai{0) + (1 - Kgal)aDM{0), 



(9) 
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Figure 4. The effect of the q and S parameters on the location 
of the critical curves (or effective Einstein radius). As seen in this 
figure comparing different (absolute-value) magnification profiles, 
significant changes in both the galaxy power law, q (top), and the 
smoothing polynomial degree, S (bottom), affect substantially the 
level of magnification since these parameters control the overall 
profile slope, but note, these only negligibly change the locations 
of the critical curves or their corresponding Einstein radius, since 
the mass enclosed by the different mass profiles is similar within 
this radius. In this example, published originally in our work on 
C10024 (Zitrin et al. 2009b), the radial and tangential critical 
curves are at ~ 10^^ and ~ 30^^, respectively. Since here we do not 
wish to constrain the (magnification) profile, but only determine 
the location of the peaks (i.e., the critical curves), it is clear that 
the choice of q and S is not fundamentally important, since for 
any such combination the curves can be formed at the right loca- 
tion, so that for our purposes these parameters can be maintained 
constant on typical values. For generating this figure the critical 
curves were originally constrained using multiple-images (Zitrin 
et al. 2009b), while in this work we show that the location of the 
critical curves can also be reasonably constrained by using only 
a proper M/L relation. 





Figure 5. Joint minimisation of the relative weight of the 
galaxies, KgaU ^-nd the M/L-dependent normalisation factor Kq, 
obtained by comparing the Einstein radii of the calibration sample 
as derived from detailed HST-based analyses, with the results of 
the automated analysis presented here. The top panel shows a 
map, whereas the bottom panel shows the 68.3%, 95.4%, 99% 
and 99.99% confidence levels {color countours), and the point of 
best fit (black circle). When these are jointly fit to the data, the 
two parameters minimised in the figure are strongly correlated, 
so that by fixing the relative galaxy weight, KgaU to the best- 
fit constant value, we are able to reduce the number of of free 
parameters in our modelling to one, namely, the M/L-dependent 
normalisation factor Kq. Explicitly, this is done by fitting a least- 
squares line to the points within /S.-}^ = 2.3 from the minimum 
X^, thus obtaining the relation between the two parameters. The 
residuals around the minimal values are then used to obtain 
the la errors. From this we obtain best fitting values (and la 
errors) of Kgai = 11.4 ± 0.6% and Kq = (51.6 ± 1-9)^. For 
more details see ^ 
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where Kgai is the relative contribution of the galaxy 
component to the deflection field. 

Both components of the deflection field are normalised 
by Kq, so that knowing its value enables us to approximate 
very well the overall deflection field. It should be stressed 
that although the degree of smoothing (S) and the index 
of the power- law (q) are the most important free parame- 
ters determining the mass profile, their effect on the Ein- 
stein radius is negligible. Based on the detailed analysis of 
^ 30 clusters (mentioned above and several more still un- 
published) , we note that the best-fitting parameters q and S 
show relatively little scatter among the different lenses. We 
can securely determine that the power-law q will be in the 
range 1 ^ ^ ^ 1.5, and the smoothing polynomial degree S 
will be in the range 4 ^ 5* ^ 24, with a sufficient resolution 
of Aq = 0.05 and AS = 2, in order to expand the full plau- 
sible profile range per cluster. More importantly, the exact 
choice of q and S does not affect the deduced Einstein radius 
size, which is determined by the inner mass enclosed within 
it and not by the mass profile which varies (see Figs. 1 and 
2 in Zitrin et al. 2009b). The crucial point to make here is 
that Einstein radii just constrain the enclosed mass, no mat- 
ter how the mass is distributed. This is seen very clearly in 
Figure [4] here, where we show that for many different combi- 
nations of the q and S parameters, the critical curves form 
at the same radius, given a reliable constraint. For exam- 
ple, in Figure [4] the critical curves were constrained using 
multiple-images, while our point here, in this work, is to 
show that the M/L ratio deduced from a calibration sample 
can be used as an alternative constraint, enabling an au- 
tomated SL analysis. In addition, it is therefore clear, that 
the choice of q and S parameter values is not fundamentally 
important here, and any combination, after it is calibrated 
for, should in principle yield the critical curves at the right 
location. 

With this in mind, throughout the analysis here we 
maintain q and S constants with ^=1.2 and iS^lO, which 
are typical values according to our many previous analyses. 
With q and S kept constants at these values, we now con- 
strain the fixed value of the weight of the galaxies relative to 
the dark matter, Kgai, and the overall normalisation factor, 
Kq. Having 10 clusters as a reference sample, and 2 param- 
eters to constrain, we can well determine their values by a 
joint minimisation, and in turn examine how these best-fit 
values reproduce the reference sample critical curves. Ex- 
plicitly, we perform a minimisation by comparing the 
Einstein radii of the calibration sample, deduced from de- 
tailed analyses based on HST observations and identification 
of multiple-images (see also with the results of the au- 
tomated procedure presented here, operated on the same 
clusters in SDSS data: 

N 

X = - <f ^^)V('^?)], (10) 

i 

where i goes from 1 to A/" = 10, for the ten calibration- 
sample clusters, and Gi is taken as 10% of the HST deduced 
values, which is a typical value for the uncertainties in SL 
modelling results. 

The results of this minimisation are seen in Figure 
[5] As can be seen, there is a strong correlation between the 
two parameters, Kgai and Kq^ which are degenerate so that 



many combinations of these can yield a good solution. This 
is a crucial point to make, since this correlation shows that 
indeed the number of free parameters in our modelling can 
be effectively reduced to one. By fixing the relative galaxy 
weight {Kgai) to its best-fit value, the model can now be 
constrained with one single parameter, namely, the M / L- 
related parameter Kq. To do this, we fit by a least-squares 
minimisation, a line to the minimum points (defined as 
lying within a Ax^=2.3 above the minimal x^)? thus ob- 
taining the linear relation between them. The fit is very 
good, — 0.97, reflecting the strong correlation, and from 
which the independent la errors are derived (i.e., by the 
residuals around the minimal x^ values). With this we ob- 
tain a best-fit (and la errors) relative galaxies weight of 
Kgai — 11.4 zb 0.6%, similar to the value expected based on 
our many previous HST analyses. 

However, one cannot expect the power-law lumpy com- 
ponent to represent only the galaxies, nor the smooth com- 
ponent to represent solely the DM, so that trying to assess 
the true physical weight of each component would be un- 
warranted at present. One can only know for certain that 
the combination of the two with Kgai as the relative weight, 
yields a good solution. It should also be mentioned, that 
we make a prior assumption on the range of sensible Kgai 
values, so that the critical curves are not too smooth nor 
too lumpy. This is done by inspecting the resulting critical 
curves by eye, so that roughly, the degree of "complexity" 
of the critical curves is similar to that seen in the afore- 
mentioned HST-based analyses of some of the calibration- 
sample clusters, and in agreement with the general expecta- 
tion for the (small) contribution of galaxies relative to the 
total mass. 

For the normalisation factor, we obtain in the minimi- 
sation a best-fit value of Kq — (51.6 =b 1.9)^^ (la er- 
rors). Accordingly, the M/L related coefficient, equals 
2.21 X lO^V^^i^"^ in units of [^r/cm^"'^; with ^=1.2], from 
which we can deduce the explicit typical M/L relation: 

M^<e)/LB = 8.7 ± 0.3 x 10"^L, 0^-'' [Mq/Lq], (11) 

where ^ = 1.2 here, is in radians, and L^ is the galaxy lu- 
minosity (in solar unit). For example, for a typical BCG 
as bright as 10^° L©, this yields an M/Lb value of ^ 
120 [Mq/Lq] within 3'', or e.g., - 180 [Mq/Lq] within 
5^^ Note that this is not the typical M/L per galaxy, but 
overall, a scaling which describes, per L©, the total pro- 
jected mass enclosed along the line-of-sight and within a 
cylinder of radius centred on a galaxy, and thus includes 
major contribution from the cluster DM halo along this line. 
Therefore, this M/L term is coupled to the modelling proce- 
dure applied and include some internal rescalings and com- 
pensation to various effects such as the difference in the 
depth between the usual SDSS and HST imaging, and the 
red-sequence membership definition (^, and are coupled to 
the LRG template and its possible minor redshift evolution 
(Benitez et al. 2009). In fact, here we do not explicitly take 
into account this evolution of red-sequence galaxies and their 
host clusters, so that the resulting M/L relation presented 
here may include a compensation to this effect, which al- 
though is expected to be minor (Benitez et al. 2009), would 
be interesting to probe when a larger calibration sample is 
available. 
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Figure 6. Einstein radius distribution (zs = 2) of :::::::10,000 SDSS clusters in different redshift bins. In each frame we specify the number 
of clusters in that redshift bin, and the median and mean radii of the distribution. The distributions are dominated, as expected, by 
galaxy and galaxy-group scale lenses, where the median and mean Einstein radii peak at zi ~ 0.1 — 0.2, although note that these rise 
again towards relatively higher redshifts of ~ 0.5, so that large Einstein radii {6e >~ 40^^) are more commonly deduced also for clusters 
at these redshifts (see also ^3.2| and |3.3| . 



With these best-fit values we analyse each of the SDSS 
reference-sample clusters sequentially in an automated way. 
The la errors on these parameters mentioned above, propa- 
gate typically into a minor ~ 5% error on the Einstein radius 
(of each individual cluster of the calibration sample), which 
may be too low in light of other uncertainties detailed in 
^ accordingly, a more realistic error (or uncertainty) level 
is estimated as mentioned therein. The resulting Einstein 
radii of the SDSS blind analysis are compared to the results 
of detailed analyses in Figure [l] where a very good correla- 
tion with a small scatter is found (R^ = 0.97, and deviation 
of less than < 5%). A more explicit example of the analy- 



sis results obtained by the different approaches is given in 
Figure [M] 



3 RESULTS, DISCUSSION, AND 
UNCERTAINTY 

The sample analysed in this work is drawn from the Hao et 
al. (2010) SDSS cluster catalog. As mentioned in their work, 
Hao et al. (2010) have developed an efficient cluster finding 
algorithm named the Gaussian Mixture Brightest Cluster 
Galaxy (GMBCG) method. The algorithm uses the Error 
Corrected Gaussian Mixture Model (ECGMM) algorithm 
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Table 1. Calibration-sample clusters (in arbitrary order). Columns are: cluster name; RA and DEC (J2000.0); photometric redshift and 
error from the Hao et al. (2010) catalog; Spectroscopic redshift; Einstein radii in arcseconds for Zs = 2 adopted from detailed HST SL 
analyses, by our automated method in SDSS imaging, and by the Jackknife minimisation, respectively; Richness, number of galaxies 
assigned to the cluster in the Hao et al. catalog; reference for the adopted 6^^^ value; other complementary references with relevant 
SL analyses. Spectroscopic redshift is assigned from the Hao et al. (2010) catalog if available, or else from the works mentioned below 
and references therein. Reference abbreviations are: Z=Zitrin, A., unpublished; 56^5=Broadhurst et al. 2005a; HSP06=Yid\V.o\di^ Seitz 
&; Panella 2006; L07/08/10=1.\uyous\iy et al. 2007/2008/2010; R09/10=K\c\id.id et al. 2009/2010, and additional references therein; 
U09=\]metsM et al. 2009; ZlO/lla/h = Zitrin et al. 2010/2011a/b; i:>7i =Donnarumma et al. 2011, and references therein; C12=Coe et 
al. 2012. Comments: We note that for RXJ2129 and A963, RIO obtained ~ 60 — 70% lower values than we obtain with the same HST 
data, although large part of the difference may be as a result of a different Einstein radius definition. While we simply measure the area 
within the explicit critical curves and from that calculate the effective radius, they quote the radial distance in which the average hi is 
equal to 1. Although both definitions are legit, they can only be expected to show similar values when the mass distribution is fairly 
round or symmetric. In addition, A2261 was previously analysed by us in Subaru imaging (U09) without any redshift or multiple-image 
information, and was recently thoroughly revised, as published in C12, based on new CLASH program imaging and combining several 
SL+WL techniques. According to the new photometric-redshift information for the multiple images found by our model (see C12), the 
Einstein radius seems to be about ~ 60% lower of that reported initially in U09, and ~ 30% lower than the value we had adopted here. 
We address the reader to these complementary works. 
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0.2233 


35.0 


36.3 
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74 
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(Hao et al. 2009) to identify the BOG plus red sequence 
feature and convolves the identified red sequence galaxies 
with a spatial smoothing kernel to measure the clustering 
strength of galaxies around BCGs. The technique was ap- 
plied to the Data Release 7 of Sloan Digital Sky Survey and 
produced a catalog of over 55,000 rich galaxy clusters in the 
redshift range 0.1 < z < 0.55. The catalog is approximately 
volume limited up to redshift z ^ OA and shows high purity 
and completeness when tested against a mock catalog, and 
when compared to other well-established SDSS cluster cata- 
logs such as MaxBCG (Koester et al. 2007; for more details 
see Hao et al. 2010). 

We go over the Hao et al. (2010) catalog, and apply 
the method described above (^ to each cluster, deriving 
its resulting Einstein radius and mass. We present here the 
results from the first 10,000 clusters analysed. In practice 
these 10,000 SDSS clusters comprise only a relatively small 
fraction (^ 20%) of the full catalog coverage, whose analysis 
results we aim to presented in a future work, once a larger 
calibration sample is available. 

3.1 Einstein Radius Distribution 

The resulting Einstein radius distribution for this sample is 
seen in Figure [g] as a function of lens redshift (for constant 
Zs =2), along with the average and median Einstein radii 
for each redshift bin, which evolve in redshift and peak at 
zi ^ 0.1 — 0.2 (for Zs = 2), as may be generally expected 
given the hierarchical growth history of clusters and the dis- 
tances involved in lensing (this is further discussed in ^3.3). 
The Hao et al. (2010) cluster catalog lists clusters with at 



least 8 members within 0.5 Mpc from the BCG. This low 
limit results in a realistic domination of galaxy-scale lenses 
{Oe of an order of a few arcseconds), which are usually not 
massive enough to form impressive lenses with large Einstein 
radii and many multiple-images. The more interesting infor- 
mation may be the higher end of the distribution at larger 
radii. The concept of the largest Einstein radius in the Uni- 
verse and the expected abundance of large lenses have been 
discussed thoroughly in the literature, and are especially of 
high interest as they teach us about the reliability of the 
standard ACDM model in predicting these extreme cases, 
as the ACDM model does not favor the formation of giant 
lenses (e.g., Broadhurst & Barkana 2008, Sadeh & Rephaeli 
2008, Zitrin et al. 2009a). 

Note that clusters with large Einstein radii are found 
also towards higher redshifts. In addition, though not in- 
cluded in this work, the largest known lens to date, MACS 
J0717. 5+3745 is at a similarly high redshift of zi ^ 0.55, 
with Oe ^ 55'' for Zs 2.5 (see Zitrin et al. 2009a). Due 
to a very shallow mass distribution in this cluster (Zitrin 
et al. 2009a), for Zs = 2 the Einstein radius will be only 
slightly lower, around Oe ^ 50'' (see also recent paper by 
Limousin et al. 2011 for new redshift information for this 
cluster). The exact number will be derived elsewhere, in the 
framework of the CLASH program. The abundance of larger 
lenses at these redshifts is caused usually (e.g., Zitrin et al. 
2011a), by a spread-out, unrelaxed matter distribution. At 
these higher redshifts many clusters are not yet relaxed and 
still undergo mergers, so that the mass distribution is al- 
ready sufficiently dense for significant lensing, but widely- 
distributed so that the Einstein radii of the different sub- 
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Figure 7. Cumulative Einstein radius distribution from ~10,000 
SDSS clusters (0.1 < z < 0.55). The cumulative distribution, and 
its upper and lower la limits, are shown in blue^ red, and green 
solid lines, respectively. Also plotted is the distribution predicted 
by the semi-analytic calculation of Oguri & Blandford (2009), 
normalised to the effective sky area of our sample (black aster- 
isks, including errors). The two distributions disagree in two main 
aspects: there is a ~ 1 — 2 orders-of-magnitude number discrep- 
ancy between them, and in addition, the two distributions have 
different slopes. The origin of the discrepancy is not clear and will 
be investigated elsewhere, although it may be as a result of dif- 
ferent mass limit, and the choice of concentration-mass relation 
and mass function used in the semi analytic calculation. Corre- 
spondingly, we find a higher abundance of large Oe clusters than 
predicted by the semi-analytic calculation. Our analysis yields 
~ 40 candidates with Oe > 40^^ (zs = 2), with a maximum of 
Oe — 69 ±12^^ (zs = 2) for the most massive cluster. Interestingly, 
this value is in agreement with the estimate by Oguri &: Bland- 
ford (2000) for the largest Einstein radius. For more details see 



structures are merged to form extended critical curves (e.g., 
Torri et al., 2004, Dalai, Holder, & Hennawi 2004). On the 
other hand, at a lower redshift, more concentrated clusters 
are those yielding larger Einstein radii, as there is more mass 



in the centre enhancing the critical area (see also J 3.3). 

The blind analysis performed here yielded initially 69 
candidates with Oe > 4:0" (zs — 2), many coincident with 
various Abell or MACS clusters. We visually inspect each 
of these clusters and find that some are boosted by various 
effects detailed below (we omit these clusters from our fur- 
ther analysis), but infer that at least about half of these are 
most likely real giant-lens candidates, with a maximum of 
6*6 ^ 69 db 12'' (zs = 2) for the most massive candidate. We 
direct the reader to works by Hennawi et al. (2007a) and 
Oguri & Blandford (2009) which have investigated in detail 
the Einstein radius abundance on various scales, based on 
simulations and ACDM expectations, and taking into ac- 
count triaxialities which induce a prominent lensing bias. 

Our realistic, observationally-based results free from 
lensing bias, are compared to some such expectations ex- 
plicitly in Figure [7| where we plot the cumulative distri- 



Figure 8. To assess the difference from the semi-analytic expec- 
tation by Oguri &; Blandford (2009; see also Figure [tJ, we com- 
pare the width of the tails for Oe > 10'\ which is the lower limit 
taken in their calculation. The histogram shows the results from 
1851 SDSS clusters following our analysis; the filled-circles curve 
shows the all-sky distribution from Oguri &; Blandford (2009), 
and the open-circles curve shows the same distribution normalised 
to the same sky area as our distribution. Both distributions are 
(semi) log-normal although with two main differences. The Oguri 
&; Blandford (2009) distribution has a width of cr = 0.1448 (in 
Log(Oe))i while our distribution shows a slower (or wider) de- 
crease, with (7 = 0.2436. Also, the overall number of clusters in 
their analysis for the same sky area, is much lower. For more 
details see also Figure |7| and ^ 



bution of clusters above each radius with the expected la 
errors, propagated from the errors on the best-fitting pa- 
rameters as described in §3.2| Note, the lower limit shifts 
the maximal Einstein radius from Oe — 69" to Oe ^ 57'' 
(zs = 2), close to that of the largest known lens, MACS 
J0717.5+3745 (Zitrin et al. 2009a). We note that Oguri k 
Blandford (2009) who examined in detail the Einstein ra- 
dius distribution based on semi-analytic expectations, have 
derived maximal Einstein radius values of ^ 60", but these 
as shown in their work are very susceptible to the cosmolog- 
ical parameters in general and to ag in particular, and can 
reach (within the 3a confidence) values that are nearly twice 
as high. Their expected distribution, scaled to the same sky 
area as our sample and with WMAP7 parameters (Komatsu 
et al. 2011), is overplotted in Figure |7| Aside for an agree- 
ment between their expected largest Einstein radius and the 
largest lenses found in our analysis, the two cumulative dis- 
tributions clearly disagree. Although normalised to the same 
effective sky area, there is a ^ 1 order-of-magnitude num- 
ber difference for small Einstein radii, which reaches a ^ 2 
orders-of-magnitude difference for higher Einstein radii, so 
that in addition, the two distributions have also different 
slopes. The origin of the discrepancy is not clear, but part of 
the difference may be due to a different (lower) mass limit 
probed by the two methods. In addition, the effect of the 
concentration-mass (c — M) relation and the chosen mass 
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function used in semi analytic calculations should clearly 
have a strong influence on the resulting distribution (e.g., 
Duffy et al. 2008, Maccio et al. 2008, Prada et al. 2011; for 
differences among various c — M relations), as higher con- 
centrations entail higher inner mass and Einstein radius. We 
leave the examination of how these may influence the cumu- 
lative distribution, for future work. 

To assess the difference from the semi-analytic expec- 
tation by Oguri & Blandford (2009; see also Figure [7|, we 
compare the width of the tails for Oe > 10' \ which is the 
lower limit taken in their calculation. As can be seen in 
Figure [s] both distributions are (semi) log-normal although 
with two main differences. The Oguri & Blandford (2009) 
distribution has a width of a = 0.1448 (in Log{Oe)), while 
our distribution shows a slower (or wider) decrease, with 
a — 0.2436. This difference is commensurate with the dif- 
ferent decline of the cumulative distribution seen in Figure 
I?! Also, the overall number of clusters in their analysis, for 
the same sky area, is much lower, although this may be, 
as mentioned, entailed by the different mass limits probed 
by each method, which will be checked in future work. As 
mentioned above, one of the most interesting aspects of our 
study, namely the largest Einstein radius, is in good agree- 
ment with the estimation by Oguri & Blandford (2009). 



3.2 Uncertainty and Error Estimation 

Various factors of error should be taken into account when 
addressing the results reported in this work, though these 
are mostly statistical and therefore due to the extensively 
large sample are not significant. In addition, these uncer- 
tainties arise mainly from the data themselves, so that ap- 
plying the method presented here to higher-end, dedicated 
cluster-survey data (e.g., the expected J-PAS survey; Moles 
et al. 2010), should produce much cleaner results. 



3.2.1 Possible Factors of Error 

The first error factor we have investigated is the lens 
photometric-redshift error. The typical photo- 2; uncertainty 
for the sample BCGs (by which we determine the lens dis- 
tance) is 0.015. A ^ 10% error in the lens redshift can be 
translated into a noticeable (> 10%) difference in the mea- 
sured Einstein radius, and only about half of the sample 
has spectroscopic redshifts for the BCG, thus the results for 
some of the clusters are affected by this error. As can be seen 
in Hao et al. (2010), the photometric redshifts were tested 
against the spectroscopic redshifts where possible, yielding 
a very tight relation strengthening the confidence in them. 
We have tested the effect of the photometric redshifts on 
the calibration sample, and regenerated Fig. [l] based on the 
photometric redshifts (instead of the spectroscopic redshifts; 
see also Table [T]). Only slight differences are seen and the 
overall scatter remains essentially the same. To assess the 
effect of the photo- 2; error more quantitatively, we analysed 



a sample of 500 random clusters (detailed in J 3.2.2 ) with the 
catalog photometric redshifts, and then repeated the analy- 
sis by photometric redshifts drawn randomly from a normal 
distribution centred on the catalog photometric redshift for 
each cluster, with a width of a = 0.015 (which is the photo-z 
error quoted in Hao et al. 2010). From this we indeed obtain 
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Figure 9. Einstein radius (log) distribution, from ~10,000 SDSS 
clusters (0.1 < 2; < 0.55). The sample has a log-normal distribu- 
tion, with {Log{9e")) 
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Figure 10. Luminosity distribution (log). Plotted is the his- 
togram of total B-band solar luminosities for each cluster, i.e., 
the sum of all cluster member luminosities. We converted the 
SDSS r-band luminosities to B-band (Vega) luminosities by the 
LRG template characterised in Bemtez et al. (2009). 



a low uncertainty of only 1.15% on the cumulative Einstein 
radius distribution, and for the differential (log-normal) dis- 
tribution, differences of only 0.38% and 0.92% on {Log{Oe)) 
and cr, respectively. 

Second, the SDSS imaging is shallower than typical 
HST imaging dedicated to SL analysis. Correspondingly, 
and supplemented by the conservative cluster-finding algo- 
rithm, some of the cluster members are overlooked and often 
not associated with the cluster, and only the brighter galax- 
ies are incorporated. Luckily, these are also the more massive 
galaxies and thus the effect on the lens model is minor. In 
addition, the inclusion of (less-massive) cluster members is 
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known to affect locally the shape of the critical curve, but 
not to change their overall size (e.g., Flores, Mailer & Pri- 
mack 2000, Meneghetti et al. 2000). The constant Kq which 
was iterated for (which includes the M/L ratio) is probably 
boosted by the relative loss of galaxy mass-representations 
in our modelling of the SDSS catalog. This, however, can 
be very well assumed to be a relatively constant ratio and 
thus, along with the M/L intrinsic scatter, does not affect 
substantially the results, as can be seen in the calibration 
sample comparison, where a clear consistency is found. This 
may not at all be surprising, since the modelling here is 
based on simple physical considerations: it has been well es- 
tablished that light approximately traces mass, and clearly, 
a reasonable M/L relation can be incorporated. 

Another factor of possible contamination is the lower 
resolution of SDSS images compared with typical lensing 
images by HST. This we find may result in local overesti- 
mation of the BCG, if another cluster member is found too 
close to the BCG core to be resolved, thus boosting the Ein- 
stein radius (and mass), especially for higher-redshift lenses. 
The reason is that the smoothing procedure, or the (2D poly- 
nomial) fit, is dominated by the BCG. Therefore, although 
the lumpy (galaxy) component in such scenarios should not 
have a substantial effect, the smooth component will be over- 
boosted in the middle (since the BCG would be too bright), 
thus pushing outwards the Einstein radius. However, this 
chance alignment is naturally not too common, and in any 
case will affect mostly the lower end of the distribution, i.e., 
clusters with small Einstein radius that is dominated fully 
by the BCG. Clusters with large Einstein radii will be less 
susceptible to such contamination as the critical curves are 
not fully dominated by the BCG and contain substantial 
contribution from other massive cluster members as well. 

A known factor of systematic error considered in related 
work on various samples of (often SDSS) clusters, is the 
miss-centring of mass with respect to the BCG (e.g., Becker 
et al. 2007, Johnston et al. 2007a, Rozo et al. 2009,2010, 
Oguri & Takada 2011). In that sense, methods which depend 
on a predefined centre may be affected from a scatter in the 
location of the BCG with respect to the true centre-of-mass, 
if the prior is constantly assumed to lie at the very centre 
of the cluster. However, in our method, there is no need to 
predefine the exact centre-of-mass. The smoothing proce- 
dure we implement has the advantage of being independent 
from a predefined centre, and the result is ultimately deter- 
mined simply and directly by the galaxy (light) distribution. 
In fact, this has enabled us to find various such shifts be- 
tween the BCG and the centre of (dark) mass (e.g., Zitrin 
et al. 2009b, Umetsu et al. 2011a). 

However, if the GMBCG catalog itself has misidentified 
a galaxy as the BCG, which we use as the centre-of- frame 
for our analysis, this may entail a shift in the analysed field, 
so that in principle some relevant galaxies may lie outside it. 
Nevertheless, since the Hao et al. (2010) catalog considers 
galaxies within 0.5 Mpc, even for the highest redshift clus- 
ters of the sample {z = 0.55), this size translates into :^ 78 
arcseconds. Since the Einstein radius is determined by the 
mass enclosed within it, and since only less than a handful 
of clusters may have such a large Einstein radius (follow- 
ing the la upper limit, see Figure |7|, this may only have a 
negligible effect over the whole sample. 

It should be noted that the results for zi > 0.43 should 
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Figure 11. Total luminosity versus Einstein radius. As is evident, 
the total luminosity itself is not an accurate indicator of the Ein- 
stein radius. Following the more realistic procedure described in 
this work is necessary in order to obtain a reliable mass model 
and consequently a reliable Einstein radius distribution. 



be more cautiously addressed, as the catalog is officially 
volume limited up to this redshift due to the luminosity 
cuts that require potential member galaxies to be brighter 
than 0.4L*, where L* is the characteristic luminosity in the 
Schechter luminosity function. Also, for higher redshifts, the 
different red-sequence criteria {r — i instead of ^ — r, see Hao 
et al. 2010) may come in play and input some more noise, 
mostly with respect to the richness level, so that overall one 
should expect fewer members assigned for zi > 0.43 clusters 
relative to clusters below this redshift. In order to test this 
effect we repeated our analysis including only clusters in the 
volume limit of 2;^ < 0.43 and verified that only negligible dif- 
ferences are seen with regard to the Einstein radius distribu- 
tion (e.g., such analysis yields a log- normal Einstein radius 
distribution with {Log{Oe")) = 0.75 and a = 0.31, similar to 
the full sample; see Figure[9|. In addition, high redshift clus- 
ters in the calibration sample also show a satisfying result 
following the same scaling relation as lower redshift clusters 
(with a scatter of up to :^ 15% with the best-fitting parame- 
ters, or up to 5% scatter with the Jackknife minimisation, 
see ^3.2.2). Also, if we exclude these from the calibration- 



sample minimisation, the best-fitting parameters differ by 
less than 3% from those obtained with the full sample. 

We have identified within the 69 initial candidates with 
Be > 40^^ {zs = 2), several clusters that were misidentified as 
higher redshift clusters (according to their observed BCG), 
though they are most likely substructures of a foreground 
more massive (and known) cluster on the same line-of-sight. 
This boosts significantly the Einstein radius, and such cases 
as mentioned were omitted from our further analysis. Due to 
the low chances of such alignments and resulting misidenti- 
fications, the effect of this on the full sample and especially 
on the lower 6e regime, is expected to be minimal. 
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Figure 12. Consistency check for our procedure. The Einstein 
mass versus the Einstein radius of the ~10,000 clusters analysed 
here. As can be expected, the Einstein radius and mass correlate 
with a square relation, whereas some scatter is seen since nat- 
urally the lenses are not strictly symmetric. Spurious detections 
did not follow the presented relation, aiding us in excluding them 
from further analysis. 



S.2.2 Quantification of Errors and Uncertainty 

In order to assess better the amount of statistical uncer- 
tainty caused by the various factors (e.g., ^3.2.1), we have 



first examined by eye a sample of 100 random clusters from 
the catalog and the critical curves generated for them by 
our automated modelling. We found only 3 clusters whose 
Einstein radius is boosted due to an unresolved galaxy near 
the BCG, and ^15 more clusters with some galaxies that 
by eye do not necessarily seem to have similar colors, so 
that if these are misidentified as cluster members, may be 
introducing some additional noise. Clearly, this designation 
is not purely objective, but still allows us to conclude, in 
addition to the other consistency checks we performed, that 
the overall noise level in our analysis is reasonable. As an 
additional complementary step, and regardless of the cali- 
bration sample, we searched for SDSS ~ 2 arcs found 
in Bayliss et al. (2011b) and examined how well our crit- 
ical curves could in principle reproduce these (giant) arcs. 
Although the location of the arcs is not used as input, our 
blind analysis automatically reproduces critical curves that 
pass through them as expected, strengthening further our 
automated approach. An example is given in Figure [15] 

To explicitly quantify the errors and level of uncertainty 
we perform two majors procedures. Firstly, we perform a 
"Jackknife" 



10) 



minimisation: on top of the x (^q- 
imisation with the full calibration sample (to find the best 
parameters for the blind analysis) , we perform the minimisa- 
tion 10 more times, each time omitting one cluster from the 
fit, and then analysing it with our automated procedure to 
examine how well the Einstein radius is estimated. We note, 
that by doing so the best-fit values for Kgai and Kq in each 
such iteration distribute around the best-fitting parameters 
when minimised by all ten clusters together, with values up 
to ^ 3a away. With this, we obtain that the Einstein radii 
for all ten clusters are estimated within ^ 17% of their refer- 



ence value (according to HST-based detailed analyses with 
multiple-images as input), while 9 clusters show a scatter of 
up to ^ 13%. These (only) represent how well each of the 
Einstein radii of the calibration sample can be reproduced 
individually. 

Therefore, secondly, we wish to examine how the (la) 
errors drawn from the full reference sample minimisation are 
propagated into the statistical results, since these should de- 
pend on other quantities such as, e.g., the number of lenses 
per (Einstein radius) bin. For that purpose, we analyse 500 
random clusters with the best-fitting parameter values, and 
then repeat the analysis marginalising over the la errors. 
For the differential, log-normal Einstein radius distribution, 
these result in differences of on {Log{Oe)), and r^T 

on cr, and upper and lower limits of ^ 18% on the cumula- 
tive Einstein radius distribution (Fig. [7|. We take these to 
represent the level of (statistical) uncertainty in our analysis. 

It should be mentioned, that although the sample anal- 
ysed here was not selected based on mass or arc abundance 
and thus is not biased in terms of lensing, the calibration 
sample used to determine the model parameters (Kq in par- 
ticular) consists of 10 well-known massive clusters, which 
might introduce a systematic error boosting the Einstein 
radii. Though low and moderate-mass lensing clusters are 
hard to model for comparison due to lack of multiple-image 
constraints, the calibration sample contains clusters with 
as few as 11 members, and as many as 142 members, thus 
spanning nearly the full richness range of the probed SDSS 
sample. The possible bias might be further looked into by 
comparing galaxy and group-scale lenses with known promi- 
nent arcs, often found in systematic surveys for gravitational 
arcs (e.g.. Sand et al. 2005, Hennawi et al. 2008, Kubo et 
al. 2010, Bayliss et al. 2011a,b, Wen, Han & Jiang 2011), 
which should also be useful for extending the calibration 
sample and examining further this effect. 

We note that due to the approach implemented here 
which does not use multiple-images as input, the profiles 
and magnifications are not well constrained for each clus- 
ter, and the only relevant measure which we refer to is the 
effective Einstein radius (and enclosed mass), as seen in Fig- 
ure |4] Naively, one could in principle derive the mass profile 
for each lens by simply assuming different fiducial source 
redshifts and calculate the enclosed mass by implementing 
their distance-redshift relation, but this would be premature 
at this stage, as though the parameters maintained constant 
here (on typical values) do not considerably affect the crit- 
ical curves shape and size, the mass profile is susceptible 
to these and thus a separate calibration is required for each 
source redshift based on the full reference sample. This how- 
ever may indeed be plausible, as we intend to probe in future 
work. 



3.2.3 Consistency checks 

To further verify that the data used here is reasonable for 
our purpose, especially the luminosity of clusters members 
to which our method is coupled, we perform a few simple 
self-consistency checks. The overall Einstein radius distri- 
bution is plotted is Figure [9] and is clearly log-normal in 
shape, with {Log{Oe'')) = 0.73t°;°3 and a = 0.316t° °°2- 
The luminosity distribution is plotted in Figure [lO] for com- 
parison, and for each cluster we explicitly compare in Fig- 
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ure |11| the total luminosity to its resulting Einstein radius, 
where it is importantly evident that the total luminosity 
itself is not an accurate enough measure of the Einstein ra- 
dius. Following a more realistic procedure as described in 
this work is necessary in order to obtain a reliable mass and 
Einstein radius distribution. More explicitly, as the mass is 
more concentrated than the light, one must choose a more 
concentrated representation for the galaxies (e.g., the power- 
law used here), which is then simply scaled by the luminos- 
ity. Similarly, we stressed that the DM is well represented 
by smoothing the galaxies mass distribution, which is more 
efficient in practice for the inner SL region, than e.g. assum- 
ing a symmetric mass distribution such as NEW which often 
does not allow to immediately uncover the multiple-images 
by the model. 

In Eigure[T2]we plot the enclosed mass versus the Ein- 
stein radius for each cluster. Although a tight relation is ex- 
pected directly from the lensing equations, this constitutes 
an important self-consistency check. We simply measure the 
Einstein radius for each cluster as the area enclosed within 
the critical curves (exploiting the magnification sign-changes 
to estimate this automatically), where the mass is mea- 
sured by summing the surface-density in the pixels which 
fall within the critical curves. The Einstein masses corre- 
late well with the Einstein radii, with a square relation as 
expected, and with a reasonable scatter since the clusters 
cannot be expected to be strictly symmetric. Explicitly, the 
of the fit is 0.985, and the mean scatter is lower than 20% 
(although there is an excess in the scatter of up to almost 
100% around 100 kpc for some individual clusters, probably 



due to the different factors of error elaborated in ^3.2.1). 



Also, in this consistency check, clusters spuriously assigned 
with large Einstein radii due to effects detailed above did 
not follow the expected relation, aiding us to exclude them 
from our further analysis. 

3.3 Correlations with cluster parameters 

Since the Einstein radius correlates with the mass interior 
to it, some dependence on the examined cluster parameters 
which are related to the observed mass, such as redshift, 
richness, and luminosity, can be expected. For example, in 
Figure [G] we showed the Einstein radius distribution in dif- 
ferent redshift bins, where it is evident that large Einstein 
radii are observed more frequently in the lower (z ^ 0.1) 
and higher (z ^ 0.5) redshifts of the sample, whereas in 
intermediate redshifts the mean Einstein radius is smaller. 

To further quantify this effect, and since Figure [5] shows 
that the Einstein-radius distribution is log-normal, we plot 
the mean (and width) of the log-normal effective Einstein 
radius distribution in different redshift bins. The result is 



seen in Figure 13 (top), where we also fit first- and second- 
order polynomials to the data. Although this tendency is 
only of the order of the log-normal distribution widths, the 
mean effective Einstein radii steadily decrease from ^ ^ 0.1 
to z ^ 0.45, and then increase again (see Figure 13). This 



tentative decrease of the mean effective Einstein radius with 
redshift may be related to cluster evolution. For example, 
lower redshift clusters, which have had more time to col- 
lapse, relax, and virialise, are expected to have more concen- 
trated mass distributions and thus be stronger lenses (e.g., 
Giocoli et al. 2011). On the other hand, the tentative in- 



crease of the mean effective Einstein radii towards z ^ 0.5 
may be related to more substructured mass distributions, 
whose critical curves for the several merging subclumps are 
merged together to a bigger critical curve (e.g., Torri et al., 
2004, Dalai, Holder, & Hennawi 2004), although it is unclear 
at present how prominent is this effect. 

It should be noted, however, that the cluster catalog is 
volume limited only up to z 0.43 where the tentative rise 
with redshift towards larger Einstein radii sets in. No weight 
should therefore currently be given to the two highest red- 
shift bins since they may be affected by different selection 
criteria applied for cluster detection above zi ^ 0.43. To 
test this, we examined plots of the cluster luminosity and 
richness versus redshift, as larger Oe clusters should have in 
principle higher luminosity and richness. A pronounced step 
is seen at z = 0.43, so that the majority of luminous clus- 
ters (> 2 X 10^^ Lq) is found at zi > 0.43, which most likely 
renders the rise in Einstein radii above zi > 0.43 a result of 
the different red-sequence criteria applied for higher- z clus- 
ters. Given this, we ignore the mean effective Einstein radii 
above z = 0.43, and concentrate on the indication of a de- 
crease from low redshift towards z = 0.43. 

One should quantify the effect of geometry on the ob- 
served evolution of the mean (log) Einstein radius with red- 
shift. The motivation is to check what is the contribution of 
pure geometrical effects, versus, say, evolutionary processes 
of the clusters. However, this can only be done qualitatively, 
since in order to know the geometrical dependence of the 
Einstein radius, one has to know the mass profiles in ad- 
vance. For example, for a point mass Oe oc (disdi/ds)^'^ , 
while for an isothermal sphere Oe oc (dis/ds). Generally, 
for a power-law projected mass distribution oc the 
angular Einstein radius scales with the distances as Oe oc 

{disdi/dsy^''d-\u we take the power law to be isothermal, 

w = 1, the angular Einstein radius (for constant Zs = 2) de- 
creases by ^ 25% from the zi ^ 0.13 bin to the zi ^ 0.45 
bin. For mass profiles steeper than isothermal, Oe decreases 
more rapidly with lens redshift, while for flatter profiles it 
increases (or shows a more complex behavior such as increas- 
ing and then decreasing). The observed monotonic decline 
between zi ^ 0.13 and zi ^ 0.45 seen in Figure [13] is of 
^ 40%, so that a steeper profile than isothermal {w c^i 1.5) 
is needed to fully explain this decline by geometrical means. 
Any possible significance of the redshift evolution of Oe for 
cluster evolution can thus only be assessed once the mass 
profile is known. 

We also compare the observed decline to semi-analytic 
calculations, in which standard concentration-mass relations 
and mass functions are incorporated. Quite independent of 
their detailed assumptions, such calculations yield increas- 
ing mean Einstein radii with lens redshift, opposite of our 
result. One such semi-analytic calculation we use for com- 
parison here (M. Redlich, private communication) is based 
on the Press & Schechter (1976) mass function. Effective 
Einstein radii are derived from assumed relaxed (i.e. non- 
merging), triaxial NEW haloes, adopting the concentration- 
mass relation from Jing & Suto (2002) and including only 
halos with M ^ 10^^ Mq. This calculation yields, similar 
to the result by Oguri & Blandford (2009), a larger fraction 
of smaller Einstein radii than our findings, and the mean 
effective (logarithmic) Einstein radius increases with red- 
shift, contrary to the decline we observe here. This however. 
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may be partly explained (like the discrepancy in Figure [7| 
by the lower-mass limit (more lower-mass clusters in higher 
redshifts would alleviate the discrepancy), or the choice of 
mass function and concentration-mass relation implemented 
therein. A more extensive calibration sample needs to be 
obtained before strict conclusions can be drawn. In prin- 
ciple, however, our results may help to pin down the true 
concentration-mass relation, in order to be compared with 
the evolution trends obtained in semi- analytic calculations, 
independent cluster evolution studies (e.g., Maughan et al. 
2008 and references therein), or related numerical simula- 
tions (e.g. Duffy et al. 2008, Maccio et al. 2008, Prada et al. 
2011). 

Finally, we note that the (logarithmic) mean effec- 
tive physical Einstein radius, i.e. the angular Einstein ra- 
dius times the angular-diameter distance to the cluster, is 
constant across the redshift range 0.1 < zi < 0.43 with 
{LogiOe [kpc])) = 1.25 ± 0.03. Tentatively, near-constant 
physical Einstein radii with redshift can be achieved with 
standard, e.g. NEW density profiles if the concentration- 
mass relation evolves steeply with redshift. Adopting a rela- 
tion of the form c oc M"°^(l + z)"^ (e.g. Duffy et al. 2008), 
we find that /3 — 3 is required to reproduce the trend seen 
in Figure [13] Current well-established concentration-mass 
relations extracted from numerical simulations find /3 ^ — 1 
(e.g.. Bullock et al. 2001, Duffy et al. 2008). We emphasize 
once more that these conclusions are tentative and prelim- 
inary, as the results of this work should be revised once a 
more significant calibration sample is available. A more pro- 
found investigation of the indicated redshift evolution and 
its comparison to numerical simulations would thus be pre- 
mature. 

Apart for the redshift dependence, we repeated the 
above procedure and examined the mean effective logarith- 
mic Einstein radii also in different total luminosity and rich- 
ness bins (namely, how many red sequence galaxies are as- 
signed to each cluster by the Hao et al. 2010 catalog). The 
evolution of {LogiOe")) and a for these is plotted in Figure 
13 {centre and bottom)^ respectively. A (mild, given the dis- 
tribution widths) trend is uncovered as a function of both 
luminosity and richness, so that on average, higher luminos- 
ity and richness clusters, show larger Einstein radii. These 
trends are quite expected, since richer clusters are naturally 
more luminous and more massive, correspondingly (for es- 
tablished mass-richness relations see, e.g., Rozo et al. 2009, 
Bauer et al. 2012). In addition, although not specifically 
shown here, for completeness we also examined both the 
richness and luminosity in the difference redshifts bins. We 
note that the richness is ^constant throughout the volume- 
limited redshift range, while {Log{Ltot [^o])) monotonically 
increases by :^ 0.4 from z ^ 0.1 to z ^ 0.45, as can be gener- 
ally expected from passive evolution of the cluster galaxies 
(although also here the trend is insignificant given the widths 
of the distribution in each bin, which is of the same size as 
the increase throughout the range, a ^ 0.4). 

We note, in addition, that with respect to the widths 
of the logarithmic distributions in different redshifts, lumi- 
nosity or richness bins, we do not uncover any prominent 
trend (which could withhold information on, say, the level 
of different population mix in the different bins). 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 




Figure 13. Evolution of the mean effective logarithmic Einstein 
radii with redshift (top), total luminosity (centre), and richness 
(bottom). The horizontal error bars mark the bin widths, and 
the vertical error bars mark the width of the distribution in the 
corresponding bin, a. In each plot we least-square fit a linear 
curve to the data (solid lines), where for the redshift plot (top) 
we also show a second-order fit (dashed line). The curve fits in the 
redshift and richness plots include only the first seven bins, due 
to incompleteness of the catalog at higher redshifts, governed by 
higher richness and luminosity clusters. In the total luminosity 
plot, we do not show the full luminosity range, since there are 
too few clusters to deduce a representative distribution for higher 
luminosity bins than those shown. See §3.3|for more details. 
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4 SUMMARY 

In this paper we presented an automated strong-lens mod- 
elling tool, which is used to efficiently derive the Einstein ra- 
dius and mass distributions of 10,000 SDSS clusters, found 
by the Hao et al. (2010) cluster- finding algorithm in DR7 
data. We adopt the well-tested approach that light over- 
all traces mass (with a more realistic representation of the 
galaxies and DM, see and normalise according to a typi- 
cal average mass-to-light relation established here, to obtain 
a reliable deflection field based on the distribution and lu- 
minosity of bright cluster members. 

This procedure, as we have shown in many previous 
SL analyses, is sufficient to derive the critical curves with 
enough accuracy to immediately identify many multiple- 
images across the lensing field, as the primary mass distribu- 
tion is initially well-constrained. Here we used a subsample 
of 10 well-studied clusters covered by both SDSS and HST 
to calibrate and test our analysis method, and showed that 
remarkably accurate determination of the Einstein radius 
can be made in an automated way, based on the light distri- 
bution of bright galaxies, and scaled by their luminosity. A 
tight correlation is seen between the Einstein radii derived 
in detailed analyses of HST data and using multiple images 
as input, and those from the "blind" survey tool presented 
here, operated on the same clusters but in SDSS data and 
without using any multiple-images information as input. 

This efficient modelling method enables us to present 
the first observationally-based representative Einstein radius 
distribution, based on a coherent unbiased sample of the 
first 10,000 clusters in the Hao et al. (2010) catalog, larger 
by a few orders of magnitude than the number of SL clus- 
ters analysed to date. For this all-sky representative sample 
the Einstein radius distribution is log-normal in shape, with 



(Log{ee")) = 0.73 



+0.02 



0.316 



+0.004 



and with higher 



abundance of large Oe clusters than predicted by ACDM, 
and with a maximum of — 69 =b 12^^ (zs — 2) for the most 
massive candidate, in agreement with semi-analytic calcula- 
tions. 

In addition to characterising the overall Einstein distri- 
bution, we also uncover various relations with cluster prop- 
erties listed in the probed catalog. For example, as may be 
expected, a clear relation is seen between the logarithmic 
Einstein radius distribution mean (for Zs = 2), and the lu- 
minosity or richness, so that richer and more luminous clus- 
ters exhibit, on average, larger Einstein radii. An especially 
intriguing trend is found with the cluster redshift. On aver- 
age, the mean effective Einstein radii steadily increase from 
z ^ 0.45 to 2; ^ 0.1. If real, and not fully accounted for by 
geometry (this requires knowledge of the mass profile, see 



^ 3.3 ), this may be a result of cluster evolution and relaxation 



processes, which make lower-2; clusters more concentrated, 
thus boosting the mass in the centre and thus the Einstein 
radius. 

Reexamining the log-normal Einstein distribution in 
physical scales rather than angular scales, we obtain best- 
fitting values of {Log{Oe [kpc])) = 1.418 ± 0.006 and a = 
0.30=b0.01, for the full sample. Subdividing this distribution, 
the mean effective Einstein radii are constant throughout 
the volume-limited range of the catalog (0.1 < zi < 0.43), 
{LogiOe [kpc])) = 1.25 ±0.03. 

The redshift trend tentatively seen in our results could 



for instance be explained by a concentration-mass relation 
that evolves more steeply in redshift than found in numerical 
simulations. If confirmed, the redshift evolution indicated 
here could help deriving an observational concentration- 
mass relation once a broader calibration sample is available. 

The results presented here are affected to some extent 
by statistical noise and uncertainty as detailed above (typ- 
ically ^ 18%; but it should be stressed that these un- 
certainties arise mostly from the data themselves, and not 
from the modelling method, which in light of higher-end 
data will produce much cleaner results. In fact, our SL al- 
gorithm could independently verify or at least probe, the 
cluster catalog and its cluster finding algorithm, by marking 
possible misidentified cluster candidates which do not follow 
the relations we obtained throughout this work. 

Such an efficient modelling method can also aid in actu- 
ally finding massive large lenses and many multiple-images, 
which in turn could be used to fine-tune the mass model 
and profile, especially when redshift information and prefer- 
ably high-resolution deep space-imaging data are available. 
Combined with complementary data such as weak-lensing, 
this will allow for an extensive examination of many other 
cluster properties, such as the mass-concentration relation, 
or a "universal" mass profile shape (e.g., the CLASH pro- 
gram. Postman et al. 2011; see also Umetsu et al. 201 la, b). 
Further analysis results of the SDSS cluster catalog of Hao 
et al. (2010) will be presented in upcoming work, where we 
will aim to include a larger reference calibration sample to 
validate further the results and uncertainties presented here. 
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Figure 14. Comparison of three clusters from our reference sample covered by both HST and SDSS. The left column shows for each 
cluster the resulting critical curves (zs = 2) of a detailed SL analysis in HST/ACS images and based on the incorporation of many 
multiple-images to fine-tune the mass model (e.g., Broadhurst et al. 2005a, Zitrin et al. 2010,2011b). The right column shows the results 
of the "blind" automated method described here, on the same clusters in SDSS data, based on the light distribution of bright cluster 
members and without any multiple-images information. As can be seen, the automatic SDSS analysis produces overall similar critical 
curves to those derived by the detailed previous analyses. More importantly, though some local difference is seen between the two 
methods, the overall Einstein radius remains the same, as seen also in Figure^ for the full reference sample. 
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Figure 15. Comparison of two SDSS clusters with giant arcs from Bayliss et al. (2011a, b). The green ellipses mark the (approximate) 
locations of ~ 2 giant arcs identified by Bayliss et al. (2011b) in deep Gemini data, and the arc marked with a red ellipse is seen in 
available HST data, and may be at a slightly different redshift. Note that without using multiple-images location as input, the critical 
curves from our "blind" analysis pass through the giant arc as expected, strengthening further our automated approach. 
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